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ABSTRACT 

Blazars are strongly variable sources that occasionally show spectacular flares visible in various 
energy bands. These flares are often, but not always, correlated. In a number of cases the peaks of 
optical flares are found to be somewhat delayed with respect to the gamma-ray peaks. One notable 
example of such a delay was found in 3C 279 by Hayashida et al. and interpreted as a result of steeper 
drop with distance of the energy density of external radiation field than of the magnetic energy density. 
In this paper we demonstrate that in general, depending on the respective energy density profile along 
the jet, such lags can have both signs and that they can take place for any ratio of these energy 
densities. We study the dependence of such lags on the ratio of these energy densities at a distance 
of a maximal energy dissipation in a jet, on their gradients, as well as on the time profile of the 
relativistic electron injection within the moving source. We show how prominent such lags can be, 
and what are their expected time scales. We suggest that studies of such lags can provide a powerful 
tool to resolve the structure of relativistic jets and their radiative environment. As an example we 
model the lag observed in 3C 279, showing that in this object the flare is produced at a distance of a 
few parsecs from the central black hole, consistent with our previous inferences based on the spectra 
and optical polarization properties. 

Subject headings: quasars: jets — radiation mechanisms: non-thermal — acceleration of particles 



1. INTRODUCTION 

Blazars are variable in all spectral bands and over a 
very broad range of time scales, from years down to days 
or sometimes even minutes (e.g., Ulrich et al. 1997; Abdo 
et al. 2010a; Ghisellini & Tavecchio 2008). Largest am- 
plitude variations are detected in the gamma-ray band 
(Abdo et al. 2010a). In luminous blazars, belonging to 
the class of Flat-Spectrum Radio Quasars (FSRQs), the 
gamma-ray light curves usually correlate with the optical 
light curves (Chatterjee et al. 2012), but with the peaks 
of short-term flares sometimes delayed with respect to 
each other (Wagner et al. 2009). In particular, gamma- 
ray flares leading the optical flares by a few days have 
been indicated in several blazars — PKS 1510-089 (Abdo 
et al. 2010b; D'Ammando et al. 2011), PKS 1502+106 
(Abdo et al. 2010c), AO 0235+164 (Agudo et al. 2011; 
Ackermann et al. 2012), and 3C 279 (Hayashida et al. 
2012). 

Such lags are very intriguing because, according to the 
most commonly accepted external-radiation-Compton 
(ERC) models, radiation in the optical and 7-ray bands is 
produced by electrons with similar energies and deeply in 
the fast cooling regime (i.e. electron cooling time scale is 
much shorter than the source light-crossing time scale). 
In such case, lags cannot result from the difference in 
cooling time scales, as in the scenarios considered by 
Sikora et al. (2001) and Sokolov & Marscher (2005). 
Hayashida et al. (2012) suggested that the optical/7- 
ray lags may result from different distance dependencies 
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of the magnetic energy density and the energy density 
of external diffuse radiation field. Here we explore this 
possibility in detail, studying the dependence of the lag 
properties on such model parameters as gradients of the 
energy densities, the electron injection time profile, and 
the ERC-to-synchrotron flux ratio at the distance of the 
maximum electron injection rate. 

Our paper is structured as follows. In Section 2, we 
describe our model of synchrotron/ERC flares that al- 
lows us to study the 7-ray/optical lags, and study how 
the resulting lag parameters depend on the details of 
the model. In Section 3, we apply our model to lumi- 
nous blazars with dense radiative environment, and to 
the particular case of 3C 279. A further discussion and 
conclusions are given in Section 4. 

2. MODELING THE LAGS 

2.1. Assumptions and equations 

Assuming that the source of a non-thermal flare prop- 
agates down the jet with a constant velocity, we can de- 
scribe its activity as a function of its distance from the 
black hole, r = ct'T, where t' is the time measured in 
the source co-moving frame and T is the bulk Lorentz 
factor of the jet. We investigate light curves produced 
by electrons deeply in the fast cooling regime by the 
synchrotron and ERC processes. They are determined 
assuming that the observed time scale of light-travel ef- 
fects caused by the finite size of a source is much shorter 
than the observed time scale of the source propagation. 
Such assumption is justified by observations suggesting 
that the AGN jet semi-opening angles are smaller than 
the Doppler angles (e.g. Pushkarev et al. (2009)). The 
ERC process is treated in the Thomson regime, using the 
following approximate formula (Moderski et al. 2003): 

[^ERC.^^HE^obs) OC [7A r 7 l't|ERc]M , (1) 
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for the electron Lorentz factor satisfying 7 2 ~ (1 + 
2 )^HE/(f 2l/ ext), where ^he is the observed frequency of 
the high-energy component, z^ ox t is the average frequency 
of the external diffuse radiation field, N 1 is the electron 
energy spectrum, |7|erc oc 7 2 Mp Xt (r) is the rate of the 
electron energy losses via the ERC process, ?4 xt (r) is the 
energy density of external radiation field as measured in 
the source comoving frame, i b s = r/{VTc) is the obser- 
vation time, V = l/[r(l — wcos^obs/c)] is the Dopplcr 
factor, and z is the source redshift. 

The synchrotron flux is approximated in an analogous 
way by 



(iobs) oc [7-^7l7lsyn]W , 



(2) 



for the electron Lorentz factors j 2 ~ 37r(l + 
z)m e cv syn /[2eB'(r)], where i/le is the observed frequency 
of the low-energy spectral component, 7|syn oc 7 2 w' B (r) 
is the rate of the electron energy losses via the syn- 
chrotron process, and u' B (r) = £?'(r) 2 /(87r) is the energy 
density of the magnetic field carried by the source. As- 
suming a power-law injection function Q 1 = K(r)-f~ p , 
where K(r) has a maximum at r = r m , one can find that 
deeply in the electron fast cooling regime 



N y = 



oc 
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where we ignore radiative losses due to the synchrotron 
self-Compton (SSC) process. 
Combination of (3), (1), and (2) gives 
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For particular spatial distributions of the energy den- 
sities u' ext (r) oc r~P E and u' B (r) oc r~@ B , we obtain that 
the ERC and synchrotron light curves are 



Fercm„ e (x) oc K(x)- 
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and 
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u Lt( r m)/u' B (r m ), x = r/r m , and A^ = 



where q m = 

In order to compute the light curves, we need to spec- 
ify the electron injection function K(x). We choose a 
modified Gaussian function 



K(x) = x ■ exp 



(x - l)(x + cr 2 - 1) 



(8) 



with a 'dispersion' a and the property that K(x) — > 
for x — > 0. 

4 Note that the factor ;r _ (p~ 2 )fe/ 4 in the latter equation results 
from the fact that in the synchrotron process 7 oc 1/ yj B(r)' oc 



2.2. Light curves 

Equations (6), (7) and (8) show that the number of free 
parameters determining the light curves is: 3 for the ERC 
case (q m , A/3 and <r); and 5 for the synchrotron case (the 
two additional parameters are Pb and p). The number of 
free parameters determining the synchrotron light curve 
can be reduced to 3 by assuming the electron injection 
function index of p = 2. Noting that electrons injected 
with such an index in the fast cooling regime produce 
spectra with the index a = p/2 = 1 {F v oc v^ a ) 1 and 
that a ~ 1 — 1.4 are typical slopes of synchrotron spectra 
of luminous blazars (Ledden & Odell 1985), we consider 
this choice to be representative, and, with the exception 
of Section 2.4, we limit our further investigations to the 
case of p = 2. 

Before plotting some examples of the light curves and 
their dependence on the model parameters, we can de- 
duce the main features of the light curves directly from 
Eqs (6) and (7). We can see that: 

1. For A/3 = 0, the ERC and synchrotron light curves 
follow the shape of the electron injection function 
and both peak at x = 1 (zero lag); 

2. From the ratio of the two light curves (oc x A @) 
one can deduce that for A/3 > the ERC peak 
precedes the synchrotron one (hereafter we call it 
the positive lag), and for A/3 < the ERC peak is 
delayed relative the synchrotron one (negative lag); 

3. For q m 3> 1, the shape of the ERC flare follows 
the shape of the electron injection function, and 
therefore has a peak at x ~ 1 (it was the case 
studied in the Appendix to Hayashida et al. 2012). 
For g m <Cl, the synchrotron flare follows the shape 
of the electron injection function and has a peak at 
z~l. 

In Figure 1, we present three plots showing the de- 
pendence of the light curves on q m , A/3 and a. We can 
see that q m determines mainly the location of the flare 
peaks relative to the point x = 1, while a and A/3 de- 
termine mainly the distance between peaks and the flare 
symmetry. 

2.3. The length and prominence of a lag 
We define the lag length as 



Axiag = X syn - X E RC • 



(9) 



(see Fig. 2) where xerc and x syn are the locations where 
the ERC and synchrotron peaks are produced. The de- 
pendence of Axi a g on q m , A/3 and o is presented in Fig. 3. 
One can see that for a = 0.6, A(3 = 4 and q m > 1, the 
lag length reaches the value of Axi ag ^0.4 and may be 
much higher for higher values of a. 

However, as one can see from Fig. 1, even for very 
large lag lengths the flares overlap significantly and are 
not strongly contrasted. To verify for which model pa- 
rameters are these lags reasonably prominent, we intro- 
duce two parameters to quantify the lag 'prominence' 
(see Fig. 2) : 5 



C: 



ERC 



^l/,He(^ERc) — -F^,HE(^syn) 



(10) 



-FiaHeOeeRc) 

5 Note that our choice of description of the lag prominence with 
two parameters is dictated by practical reasons. Because of typ- 
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and 

n f^LEpgsyn) ~ -F^LE (^ERC ) , n 

L'syn — ^ 7 \ • (,11 J 

-fV.LE^I^synJ 

The dependence of these indicators on the model param- 
eters is presented in Fig. 4. We find that both measures 
of the lag prominence increase with increasing A/3 and 
increasing a. The dependence on q m is more complex 
for q m < 1 the lag prominence decreases with increasing 
q m , while for q m 3> 1 the lag prominence depends on g m 
very weakly. 

2.4. The case of p ^ 2 

In Figure 5, we show the dependence of the lag length 
on the injected electron energy distribution index p. We 
can see that regardless of the value of A/3, the lag length 
slightly decreases with increasing p. This means that the 
absolute value of the lag decreases for positive lags and 
increases for negative lags. This can be easily understood 
from the existence of a factor ;e~(p~ 2 )' 3b / 4 in Eq. (7). 

3. APPLICATIONS 

We have studied basic features of lags between the 
ERC and synchrotron light curves produced in the fast 
cooling regime. Our model shows that the presence of 
such lags is made possible by a non-zero value of A/3 
which corresponds to a significant change of the Comp- 
ton dominance parameter q = u' c ^ t /u' B cx r~ A/3 within 
the distance range of enhanced activity of the propagat- 
ing source. In a quasar radiative environment the change 
of q with distance is expected to result mainly from 
the strong stratification of the external radiation field 
sources, being connected with broad-line region (BLR) 
and hot-dust region (HDR). Changes of q with distance 
can be deduced from Fig. 6, where radial distributions of 
respective energy densities are schematically presented 
(see also Sikora et al. 2009). In this figure we marked re- 
gions where q decreases with a distance (region "I" ) and 
regions where q increases with a distance (region "II"). 
Hence, positive lags (t 7 < i op t) are expected to be pro- 
duced in regions "I," and negative lags are expected to 
be produced in regions "II." 

Since optimal conditions for production of prominent 
lags are such that the distance range of enhanced source 
activity determined by x ~ a coincides with the region 
of largest changes of q, one may expect that observed 
lags, A£i a g, are produced at distances r m ~ Ari ag ~ 
cAtiagXT/ (1 + z). Hence, when observing a lag of Ati ag , 
one can estimate r m and locate the source activity rel- 
ative to tblr and thdr, which values can be estimated 
from the broad line luminosities and the hot dust sub- 
limation temperatures (see Sikora et al. 2009 and refs. 
therein). Unfortunately, due to limited effective area of 
present 7-ray instruments (such as Fermi-LAT) it is pos- 
sible to construct light curves with data quality only suf- 
ficient to verify the presence of correlation between op- 
tical and 7-ray data and estimate lag length only with 
time resolution longer than a day (with exception of the 
very brightest 7-ray sources). Noting that in luminous 

ically poor time sampling and large errors in high energy light 
curves, such definition enables us to estimate the lag contrast in 
a more precise way than using a single parameter, e.g. the flux at 
the intersection of two normalized flares. 



blazars tblr ~ 10 18 cm, and that the observed time of 
the propagation of the source over the distance range of 
this order is ~ 1.7(1 + z)/(r/20) 2 days, presently only 
lags produced in the HDR have a chance to be probed 
observationally. 

To verify the applicability of our lag model, we use it to 
match the data on blazar 3C 279 presented by Hayashida 
et al. (2012) (see also Abdo et al. 2010d). A light curve 
analysis presented in their paper (R-band vs. 7-rays 
of energies > 200 Me V) with the discrete-correlation- 
function method showed a peak at +10 days, indicating 
a presence of a positive lag of that length. We used their 
data from the time interval MJD 54870-900 including 
representative strong flares in optical/7-ray bands, and 
made an attempt to model them with the flare profiles 
given in Sec. 2.1. The data and the model light curves 
are presented in Fig. 7. The parameters of our model 
are: q m = 30, j3 B — 2, (3 E = 5, a = 1, r m = 4.4 pc 
and r = 10. These results are consistent with the ERC 
model parameters used by Hayashida et al. (2012) to fit 
the broad band spectra of 3C 279. We would like to em- 
phasize that the presence of a +10 days lag in the case 
of 3C 279 is just a possibility as the used analysis might 
have been misleaded by poor data quality. Our choice 
of the strong flares from the whole lightcurve presented 
in Hayashida et al. (2012) is therefore just an example 
indicating that the model can be applied to a real data 
with very reasonable parameters. 

4. DISCUSSION AND CONCLUSIONS 

In rarefied magnetized plasmas of astrophysical non- 
thermal sources, dissipative processes often generate 
short-term bursts of relativistic electrons, which in turn 
produce synchrotron and inverse-Compton flares. Pro- 
vided that the energy density of seed photons, measured 
in the source co-moving frame, is dominated by external 
sources, the inverse-Compton radiation is dominated by 
the ERC process. For quasi-uniform energy densities of 
magnetic and seed radiation fields across the source, and 
for electrons energies such that the time scale of their 
radiative cooling is much shorter than the time scale of 
the burst, the synchrotron and ERC flare profiles will 
follow the electron injection function shape almost ex- 
actly, without any lags. At the same time, the flux ratio 
of the two spectral components will approximately be 
scaled by the ratio of energy densities of seed radiation 
field to magnetic field, as long as the inverse Compton 
scattering proceeds in the Thomson regime. 

However, in blazars, the outbursts are produced by 
sources (e.g., internal shocks, magnetic reconnection do- 
mains) which propagate down the jet with relativistic 
speeds. As we illustrated in Fig. 6, in FSRQs, the en- 
ergy density distribution of external radiation field is sig- 
nificantly stratified, while the magnetic field intensity is 
expected to monotonically drop with a distance due to 
the jet divergence. As our equations in Section 2 show, 
in such an environment, for sources propagating down 
the jet, optical light curves produced by the synchrotron 
mechanism differ from 7-ray light curves produced by 
the ERC process. In particular, if the source is burst- 
ing within the regions marked on Fig. 6 by "I" or "II," 
we will observe flares which are delayed with respect to 
each other, with the 7-ray flares preceding the optical 
ones if the burst takes place at tblr < Arburst <C thdr 
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and Arburst ~3> ^hdr, and with the optical flares pre- 
ceding the 7-ray ones for Ar burst < t B lr and t B lr < 

Afburst < fHDR- 

Our scheme does not include the presence of a warm 
dust, which together with the hot phase, produces spec- 
tra extending up to 20 /xm - 30 fim (see, e.g., Fig. 7 
in Shang et al. 2011). The observed spectra indicate 
a distribution of dust extending to much larger distances 
than thdr- In this case, the decrease of the external 
radiation energy density does not necessarily have to be 
faster than the decrease of the magnetic energy density. 
However, frequently observed excesses of the NIR radi- 
ation over the MIR component extrapolated to shorter 
wavelengths (Mor et al. 2009; Leipski et al. 2010) suggest 
that there can be two separate thermal components from 
dust, the hot one produced by graphite and located in 
the outer portions of the BLR, and the cooler one pro- 
duced by both graphite and silicate grains and enclosed 
within clumps of the molecular tori extending up to tens 
of parsecs (Mor & Netzer 2012; Kishimoto et al. 2012; 
Roseboom et al. 2012). Hence, a confirmation of pre- 
dominantly positive optical/7-ray lags of a few days can 
independently support such a stratification of the dust 
distribution. 

Obviously, the model proposed and studied by us, 
aimed to explain the observed lags between optical and 
7-ray flare peaks is not a unique one, but is the only 
one which involves single dissipative events. This model 
can be verified in future against two-dissipative-zone sce- 
narios, e.g. those involving reconfinement plus reflection 
shock structures, or two internal shocks formed indepen- 
dently at different locations (Katarzynski & Ghisellini 
2007), by future more detailed and better resolved light 
curves. 



• Future more detailed and better resolved light 
curves should allow to discriminate whether the op- 
tical - 7-ray lags can be explained by single dissi- 
pative events, or whether two spatially separated 
dissipation zones must be involved. 
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We studied the dependence of synchrotron-ERC lags 
on model parameters which determine the distribution 
of magnetic and external radiation fields, their energy 
density ratio at the time of the electron injection peak, 
the dispersion of the injection function, and the slope of 
the injected electron spectrum. Our main results can be 
summarized as follows: 

• Variations of q = u' cxt /u' B with a distance result in 
different ERC and synchrotron light curves; 

• Electron injection bursts produced within the re- 
gion of the monotonic changes of q lead to the pro- 
duction of the lagged flares, positive (7-ray first) 
for dq/dr < 0, negative for dq/dr > 0; 

• The 'prominence' of such lags increases with in- 
creasing burst dispersion <j and increasing gradi- 
ents of q. The most prominent lags are expected 
when the burst length coincides with the region of 
largest gradient of q; 

• Limited time resolution of 7-ray data currently al- 
lows only searches for lags larger than a few days, 
which corresponds to distances of the order of, 
or larger than, the sublimation radius of the hot 
dust. As an example, we showed that our model 
light curves can be matched to the observed 7- 
ray/optical flare of 3C 279 in February 2009, when 
a 10-day lag was seen (Hayashida et al. 2012). 
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Figure 6. Schematic view of distance dependence of magnetic and external radiation field energy densities as measured in the source co- 
moving frame (see Fig. 3 in Sikora et al. 2009). Label T marks the regions where positive lags are expected, whereas label 'II' corresponds 
to regions where negative lags are expected. 




Figure 7. The February 2009 7-ray/optical outburst of blazar 3C 279 (from Abdo et al. 2010d; Hayashida et al. 2012) with a time lag of 
£b 10 days in normalized units. Lines show our analytical model (dotted line — ERC flare; solid line — synchrotron flare) for parameter 
values q m = 30, /3_b = 2, Pe = 5, cr = 1, r m = 4.4 pc and Y = 10. 



